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Where  Solitons  May  Be  Hiding  in  DNA  and  Their  Possible  Significance 

in  RNA  Transcription 

E.  W.  Prohofsky 

Department  of  Physics 
Purdue  University 
West  Lafayette,  IN  47907 


We  find  that  the  hydrogen  bond 
stretch  bands  of  the  double  helix  ap¬ 
pear  to  be  nonlinear  enough  to  support 
solitary  wave  energy  concentration. 
Coupling  this  fact  to  predictions  of 
our  selfconsistent  theory  of  helix 
melting  gives  rise  to  speculations  of 
a  mechanism  for  base  pair  melting  in 
KNA  transcription  which  is  consistent 
with  the  known  energy  needs  of  tli.it 


process. 
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Several  recent  papers  have  discussed  the  possibility  of  solitary  waves  in 
double  helical  DNA  associated  with  the  helical  acoustic  modes  (1).  The  work 
was  stimulated  by  observations  in  DNA  of  resonant  microwave  absorption  (2). 

Papers  have  also  appeared  which  have  disputed  the  theoretical  predictions  of 
the  likelihood  of  such  solitary  waves  (3)  and  papers  have  also  appeared  which 
challenge  the  experimental  observations  (4).  Solitary  wave  phenomena  require 
nonlinear  effects  and  both  experimental  observations  (5)  and  theoretical  analysis 
(6)  indicate  that  helix  acoustic  modes  are  quite  linear  in  their  dispersion  and 
amplitude  dependence.  In  our  theoretical  studies  of  helix  lattice  dynamics  of 
DNA  we  have  found  that  a  different  region  of  the  DNA  spectrum  is  fur  mote 
nonlinear  (7)  and  seems  to  satisfy  the  quantitative  requirements  for  solitary 
wave  effects. 

The  highly  nonlinear  modes  are  the  hydrogen  bond  stretch  modes  at  It  - 
120  cm  These  modes  have  been  observed  in  both  low  frequency  Raman  x  mo 
ing  (8)  and  Fourier  transform  infrared  absorption  (9)  experiments.  In  b.  t h 
rases  the  modes  have  been  seen  to  soften  in  frequency  with  in<  teas  tin-  t «  ■■■  •  . 

true.  We  have  predicted  theoretically  that  these  modes  are  temper  ituti 
in  agreement  with  observations  (10). 

The  same  vibrational  modes  have  played  a  central  tele  in  the  tin  ,  t 
helix  melting  that  we  have  developed  based  on  helical  latt  i<e  dvnami  *■  .  n 


I  i  mge  i  at  lire  dependence  discussed  above1,  inupled  with  lint" 
I.  id  naturally  to  hydrogen  bond  breakup  and  helix  ,  .  1  t  i  no  . 
pel  v(iK')  .polyfdt'l  p  i  ed  t  c  t  me  I  t  i  ng  at  a  I  enipe  i  a  l  u  r  e  in  lair 
el  i  •  .  t  v  a  t  i  mis  .  In  tbit  work  all  the  pa  rul.ie  l  e  t  S  we  t  e  t  i  t  t  ed 
i  mi  i  en:pi  •  i  a  t  ure  and  no  parameters  were  titter!  to  ratting. 
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Helix  melting  is  an  essential  element  in  many  biological  processes.  For 
example,  it  is  essential  in  the  process  of  RNA  transcription.  Recent  reviews 
of  this  process  reaffirm  the  belief  that  energy  must  flow  into  the  RNA  polymerase 
complex  from  its  surroundings  (12).  The  ability  of  solitary  wave  effects  to 
concentrate  excitations  in  these  bands  would  help  in  understanding  the  problems 
witli  energy  flow  in  these  processes.  This  biological  role  for  nonlinear  effects 
is  in  line  with  the  original  suggestions  of  solitary  wave  significance  in  bio¬ 
logical  processes  (13). 

In  this  paper  we  discuss  our  helical  lattice  dynamics  results  for  the 
mi, linear  modes  and  show  how  these  results  can  be  used  to  extract  parameters 
Hided  to  evaluate  the  possibilities  of  nonlinear  concentration.  We  then  use 
r  In'  ilieorv  oi  melting  developed  to  give  order  of  magnitude  estimates  of  the 
l  I  i t  mles  needed  to  induce  base  pair  melting,  and  discuss  the  significance  of 
1  ,  ■ e  results  ter  KNA  transcription. 
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!  Stretch  Modes  in  the  Double  Helix 


Is  ’Mil  st  ml  d  the  vibrational  modes  of  the  double  helix  we  have  found 

i  1 1  I  sills  wiiih  ■  .in  be  char  act  et  i  zed  as  vibrations  that  compress  the 

I-  "'is  si  t  u  i  is)  ha  ;e  pairs.  There  are  overlapping  bands  for  all  combi  no - 
■  1  e.i  |  mis  os I  seveial  p.e;s  hands  in  this  region  should  exist  tegard- 

:  i  -i  .i  i|  u  a.  e  tlu).  Dispersion  curves,  i.e.,  the  change  in  frequency 
•  ;  a  .  iii  shewn  tor  several  such  bands  in  Figure  1  (15).  Theta 
I  •  i  !  I  el  the  d  I  sp  I  .  e  else  n  t  s  f  i  on  one  unit  cell  to  the  next  .  'I  he  so 

1  1  '  s-  !'i  i  ’  'I  i  n  i  I'm  i’.:  t  i  eijiieiii  i  es .  Tn  one  i  use  the  Imal  minimum  is 

ildd  I  i  i  oi,  in  ether  euros  its  at  a  /one  edge.  I  he  i  ill  ves  could 
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be  extended  about  either  end  i.e.  the  total  6  range  can  run  from  -180°  < 

0  •  180°  or  0  <  0  <  360°.  The  curves  would  have  symmetry  about  the 
extension  as  the  positive  and  negative  phase  shifts  are  connected  by  time 
reversal  invariance.  Correlating  frequency  with  energy, and  phase  shift  with 
wnvevector  and  hence  momentum  for  these  excitations, one  can  relate  these 
curves  to  particles  with  energy  vs.  momentum  relations  approximated  by 


The  effective  miss  m  is  then 


if  a  is  the  lattice  spacing 
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where  0  is  in  radians  and  Aw^  is  the  change  in  frequency  due  to  dispersion 
ci Toi  ts  for  a  shift  in  tlio  phase  f  r<>m  cell  to  cell  of  Lf  .  In  the  solitary 
wave  .analysis  by  Davydov  (13)  the  wave  equation  for  such  a  system  is  given 
by  the  ni'ii  I  i  near  Si  breed  i  nger  Ffpi.it  ion  in  one  dimension 
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this  becomes 
^  3x 
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which  has  the  solitary  wave  solution 
i  (2kx  -  cn  ) 
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where  V  is  the  velocity  of  the  centroid  of  the  solitary  wave.  Concentration 
of  the  excitation  occurs  in  the  square  wave  excitation  case  (13)  for  g  >  g 
where  „ 


cr 


71 

u ~i 
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and  •  is  the  extent  of  the  initial  excit.it  ion  along  x. 

The  t  err.  (i\  ,  J  “  is  a  measure  of  how  the  energy  is  lowered  by  increased  con- 
.  nitration  of  the  probability  amplitude.  in  our  sel i consistent  temperature 
dependent  lattice  calculations  we  have  found  that  the*  hvdrogen  bond  stictch 
ii.-dee  ate  . .  i  d  1. 1 11:101 1  i  i '  ,  i.e.  t  liev  1 1  v.e  l  in  f  roqiiont  v  as  the  excitation  1  i  ve  1 
a  tempi  ratine  im  tease,  l-.o  have  I  mind  this  sliift  in  frequency  to  he  in  t'nueh 
leieement  with  the  observed  !  i  eqiieiii  \  shift  from  line  temperature  to  room 
'  i  talari  (10).  hi  ill-  I  i  1  '  i '  in  e'ettV  with  ini  I  '  i‘-i'  in  I  HOI  I'l't  rat  i  CO' 

|  1  i  t  udi  j  i  i >v  i ■  i-  s  t  hi  si  1  1  t  i  | -p  i  it /  i ■  ’  i  i  .  t  i.h  i  i  h  i  d i  in  i  s  .  in  i  n i  i  a t  :  on  o )  t  l  , 

i  ■  i  i  t . '  t  i  i  a'. .  In  I  . .  '  ’  '•  o  .>w  a  i  1 1 1  I  it  ion  f  I '» >  ot  i  in-  :  i  i  :  .a  a  v  ■  ■  1  .  •  r  1 1 

■ 1 ,  ■ .  i  i  i  |  -’ll  I  de  .  *  |  ■  i  *  ]  \  (  i l  )  •  •  a.  I  I  ai  .  ■ :  t  .  1 1 ■ .  i  1 1  .  ■  i  ‘  i  •  i  t  .  ,  .  |  .it  ’  ,i 
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room  temperature.  For  vibrational  systems  the  mean  wave  amplitude  squared 
is  proportional  to  the  temperature.  Tie  frequency  is  seen  to  drop  with 
increasing  temperature  and  hence  mean  squared  amplitude. 
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We  can  estimate  C,  for  particular  excitation  modes  from  Figure  2.  In 


Kquation  3,  |  if  ]  is  a  normalized  nondimensional  squared  wave  amplitude.  We 

i  2  AT 

associate  A  |  [  with  —  for  regions  of  curve  in  Figure  2.  We  also  associate 

the  appropriate  AE  with  fiAw^  over  that  AT.  is  the  change  in  frequency 

due  to  anharmonici ty  not  dispersion.  Then: 

Am 


G  =  _ A_ 

fi  AT/T 


(9) 


I'sing  this  and  Eqs.  (2)  in  Eqn.  (4) 
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ma 

g  =  --- x  - 


, . .  N  2  Aui 
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Fur  the  85  cm  band  from  Fig, tires  i  and  2  at  1  r  room  temperature  we  estimate 

r  t  10.  For  reasons  discussed  in  the  next  section  we  estimate  P  ,  the  minimum 

run 

initial  size  of  the  excitation  to  be  2.5  in  units  of  a.  For  this  value 

g,  ~  1.  One  can  expect  nonlinear  concentration  of  excitations  in  these  hands 
<  ■  r 
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UNA  l.vbi  id  1. 


is  split,  a  I  IN’ A  base  pa  i  i  ti  .  s  I  t  i  am  j  s 


a  polymerization  event  in  the  RNA  backbone  takes  place.  As  a  result  of  the 
backbone  chemistry  n  pyrosphophate  group  is  removed  from  the  nucleoside 
triphosphate  (NTP )  which  is  the  source  of  the  added  RNA  section.  The  pyrophos¬ 
phate  group  is  later  hvdrolized  but  not  at  a  site  directly  connected  to  the 
other  operations.  The  overall  energy  flow  in  the  transcription  complex  is 
such  that  on  the  average  0.5  kcal/mole  transcribed  must  be  absorbed  by  the 
complex  from  its  surroundings. 

For  the  purposes  of  this  discussion  we  will  explore  ways  in  which  energv 
introduced  into  the  complex  can  bring  about  base  pair  melting.  The  free  energv 
needed  to  melt  a  (1C  pair  is  generally  accepted  to  be  3.5  kcal/mole  and  that 
lor  an  AT  pair  1  kcal/mole  (12).  If  inflow  of  this  amount  of  energy  occurred 
the  net  energy  requirements  of  transcription  would  easily  be  met.  The  reason 
to  consider  this  form  of  energy  transfer  to  the  transcription  complex  is  that 
we  believe  it  would  involve  the  nonlinear  hydrogen  bond  stretch  modes.  In  our 
theory  of  melting,  we  calculate  the  vibrational  fluctuations  I)  on  each  hydrogen 
bond  where  I)  =  •  u u  •  and  u  a  is  the  displacements  of  atoms  f rom  equilibrium 

due  to  phonon  amplitude.  The  se 1 1  eons i stent  softening  of  the  bond  and  thermal 
expansion  are  calculated  based  on  D.  The  softening  and  expansion  lead  to  a 
'•'ell  i  up,  oi  the  hydrogen  bonds.  In  tiiis  sense  D,  which  grows  with  excit.it  ion 
1  •  '■  ■  1  ,  drives  the  melting,.  In  our  model  where  we  looked  at  melting  adjacent  to 


i 
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m  I  t  ed  I  ise  paii  in  pn  I  y  (d(  )  •  po  1  y  ( ciC )  we  found  that  when  a  hydrogen  !  ■  :id  ins 

0.0,0  I". ■  1  I  i  in’,  i n  .  u r e.l  (  !  !  )  . 

■  *  l.o  i  n  vi  r,t  i  :  t  ed  t  he  .  out  t  ibut  ion  ol  individual  vibrational  ha -an.  1 I  ( i.  1 

t  in  !  <  t  a  I  I  1  is  t  ua  (  1 1  ai  1)  ,  i  •-  ceu  in  Figure  3 .  I  lie  h  v.l  i  o.  eii  bond  s  l  i  ,  •  l  .  ’  i  Fund 
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at  95  -  100  cm  at  340K  has  D(m)  z  5.5  x  10  A  from  the  entire  band.  Ignoring 
the  contributions  of  other  bands  one  would  only  have  to  increase  the  excitation 
of  this  band  by  a  factor  of  7.27  to  bring  on  melting  from  the  energy  in  this 
band  alone.  Since  the  thermal  occupation  is  about  2  phonons  at  this  temperature 
one  would  require  ~ 14—15  phonons  to  bring  about  melting.  The  energy  of  t  14 
phonons  is  approximately  the  3.5  kcal/mole  generally  associated  with  the  binding 
energv  of  CC  base  pairs.  The  localization  to  the  hydrogen  bond  of  these  excita¬ 
tions  is  in  the  localized  amplitudes  calculated  and  plotted  in  Figure  3. 

If  we  consider  melting  a  CC  base  pair  in  terms  of  the  needed  energy  in 
one  mode  while  all  others  are  in  thermal  equilibrium  at  room  temperature  a 

2 

dilferent  estimation  is  appropriate.  One  still  needs  a  total  D  ;  0.004  A 

2 

The  thermal  activation  of  all  modes  at  room  temperature  is  D  ;  0.028  A  (7). 

K 

2  _i 

The  needed  additional  fluctuation  is  AD  =  0.012  A  .  For  D(w)  of  the  95-100  cm 
band  I) (u0  =  0.0055  to  make  up  this  difference  the  excitation  level  would  have 
t  •  •  3.8  times  that  at  343K.  This  requires  ;  6  phonons  in  a  picket  arriving  at 
tin  base  pair  .adjacent  to  tin1  melted  region.  The  thermal  activation  level  is 


mean  1  at} 


rge  fluctuation  value  front  Bose-F inst e in  statistics  is  :4-5. 


Su.-h  a  eonienli.it  ion  of  energy  may  be  achieved  by  capturing  energy  from  the 
:  .u i  i . ’ued i ups  as  is  apparently  needed  to  explain  the  energetics  of  DNA  li.mscrip- 


I  O’e.  Sue 


h  an  i  m  it  it  ion  could  travel  in  unmelted  DNA  without  disrupting,  the 


n;;.\  .s  i  he  |>(  .  )  t  <>r  all  hands  is  less  for  base  pairs  t.ir  I  i  om  a  melted  sect  ion. 

:H  i  i  ei  lie  H)0  (m  hand  is  about  half  the  value  at  a  cel  If'')  as  for  cell  (»  1  ) 

•  |  t,i  |  In-  r.i  !  t  ,  as  . .  in  figure  3.  When  l  hi'  one  r  g,\  pi  npugal  e:  into  the  ceil 

.  i.  ,  i  <ni  i  i  |  hi  i  e  I  t  h.,n}»e  i  :t  1  oca  1  i  zed  mode  behavior  greatly  i  ne  tea  m  s  f  l;u- 

,  >  i. .  I  i  .  i ;  ,  ■  j  i  i  t  i  ■  t .  •  .mil  d  i  si  up  i  .  'he  ha.ee  pair. 


It  is  an  interesting  question  as  to  whether  the  RNA  polymerase  is  a  proper 
thermodynamic  engine  which  runs  on  a  directed  flow  of  energy  or 


is  driven  hv  thermal  fluctuations  as 


many  believe 


(Id).  If  solitary  wave  effects  can  be  shown  to  in  fact  exist  in 

I)NA  one  will  have  cornea  long  way  to  showing  that  the  RNA  polymerase  can  be  a 
proper  thermodynamic  engine.  Energy  from  the  hydrolysis  of  pyrophosphate  could 
be  absorbed  before  the  thermalizat ion  of  this  energy.  The  entire  helix  could 
act  as  an  antenna  for  the  absorption  greatly  increasing  the  fraction  of  energy 
absorbed  before  thermal izat ion .  The  processing  and  concentration  of  energy 
would  occur  as  the  energy  propagated  on  the  one  dimensional  helix  to  the  region 
of  transcription. 


In  the  earlier  section  we  estimated  8  the  initial  size  of  the  excitation 
as  being  at  least  ;  2.5  base  units  in  length.  This  is  based  on  the  idea  of  a 
planar  longitudinal  wave  or  shock  wave  from  pyrophosphate  hydrolysis  impinging 


on  a  sec 


lion  of  UNA  helix  from  the  side.  A  helical  structure  is  just  what's 


n .rled  to  convert  such  an  external  excitation  into  a  short  wavelength  excitation 


in  (lie  internal  coordinates  of  the  DNA  because  of  k-c onservat ion  selection  rules 
|'< * r  some  forms  of  excitation  a  10  fold  helix  would  select  excitations  of  10 


i  s(  pair  w  avr  I  e 


ngth  or  n  =  36  on  dispersion  curves  such  as  Figure  1.  To 


pen*  tale  a  hydrogen  bond  stretch  wave  where  compression  of  either  strand  leads 
l  i >  ,  oviprc  sion  of  i  lie  hydrogen  bonds  the  wavelength  would  be  5  base  pairs  or 

We  have  e. seriated  half  (bis  wavelength  with  the  minimum  extend  of  (he 


1  I .  I  I  r  V  \  e  r  :•  I 


i  tat  ion  ami  sel  f:  ~  2.5.  Short  wavelengths  can  leid  l< 


i  e  I...  .i  1  i  ■ .  .  1  v  avr  H  1  e  t  s  a:ul  local  i/at  ion  or  concent  i  at  ion  is  needed  to 
t  i-  i<  ui  !  v  i  elt  tie. i  pairs  ad  ••cent  the  already  melted  area.  It  is  int  ere*t  i  ng 


to  speculate  on  hew  the  existence  of  a  small  integer-fold  helical  structure 
for  DNA  interacts  with  the  need  for  energy  transfer  in  process  such  as  transcrip 
tion  in  early  evolution  of  the  system. 

Better  nonlinear  models  of  wave  propagation  more  appropriate  to  the  case 
of  dispersive  optical  modes  would  allow  a  more  quantitative  estimate  of  many 
of  the  energies  which  have  been  crudely  estimated  in  this  phenomenological  dis¬ 
cussion.  The  role  of  DNA  excitations  in  biological  processes  may  be  important 
and  -eems  likely  to  be, based  on  these  sir, pie  estimates. 

This  work  supported  in  part  by  ONR  Grant  XOOO14-86-K-0252  and  NIH  Grant 
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Figure  Captions 


Figure  1.  Dispersion  relations  for  several  optical  bands  in  poly(dG) .poly(dC) . 

The  frequency  is  in  cm  ^  and  the  phase  is  the  relative  phase  between 
neighbor  unit  cells.  In  terms  of  a  wavevector  q,  this  phase  is  qa 
where  a  is  the  lattice  constant. 

Figure  2.  Calculated  change  in  frequency  of  selected  zone  center  modes  of 
poly (dG) • poly (dC)  as  a  function  of  temperature.  For  details  see 
reference  ( 14 ) . 

Figure  3.  The  contribution  to  the  fluctuation  amplitude  by  individual 

vibrational  bands  for  the  hydrogen  bond  adjacent  the  major  groove 

in  polv(dG) -poly(dC) .  The  three  divisions  are  for  the  two  base 

pairs  adjoining  a  melted  base  pair  on  either  side  and  for  a  base 

-1 

pair  far  away  from  the  melted  base  pair.  For  the  95-100  cm  band 
the  highest  division  is  for  the  (+1)  base  pair  i.e.  3'  -*•  5'  direction 
in  the  G  chain.  The  next  highest  is  for  the  (-1)  base  pair  i.e. 
on  the  other  side  of  the  melt.  The  lowest  division  is  for  the 
(«)  base  pair  i.e.  far  away  from  the  melt.  See  reference  11  for 
details. 


f  req  ( cm 


